NEVIS-1529 



A High Statistics Search for v^v^) — > u e (u e ) Oscillations in the Small Mixing Angle 

Regime 

A. Romosan, 2 C. G. Arroyo, 2 L. de Barbara, 5 P. dc Barbara, 7 A. O. Bazarko, 2 R. H. Bernstein, 3 A. Bodek, 7 
T. Bolton, 4 H. Budd, 7 J. Conrad, 2 R. B. Drucker, 6 D. A. Harris, 7 R. A. Johnson, 1 J. H. Kim, 2 B. J. King, 2 
T. Kinnel, 8 M. J. Lamm, 3 W. C. Lefmann, 2 W. Marsh, 3 K. S. McFarland, 3 C. McNulty, 2 S. R. Mishra, 2 D. Naples, 4 
P. Z. Quintas, 2 W. K. Sakumoto, 7 H. Schellman, 5 F. J. Sciulli, 2 W. G. Seligman, 2 M. H. Shaevitz, 2 W. H. Smith, 8 
P. Spentzouris, 2 E. G. Stern, 2 M. Vakili, 1 U. K. Yang, 7 and J. Yu 3 
University of Cincinnati, Cincinnati, OH 45221 
2 Columbia University, New York, NY 10027 
3 Fermi National Accelerator Laboratory, Batavia, IL 60510 
4 Kansas State University, Manhattan, KS 66506 
5 Northwestern University, Evanston, IL 60208 
6 University of Oregon, Eugene, OR 97403 

7 University of Rochester, Rochester, NY 14627 

8 University of Wisconsin, Madison, WI 53706 

(February 7, 2008) 

Limits on — > v e (p e ) oscillations based on a statistical separation of v e N charged current 
interactions in the CCFR detector at Fermilab are presented. v e interactions are identified by the 
difference in the longitudinal shower energy deposition pattern of v s N — » eX versus v^N — > v^X 
interactions. Neutrino energies range from 30 to 600 GeV with a mean of 140 GeV, and flight 
lengths vary from 0.9 km to 1.4 km. The lowest 90% confidence upper limit in sin 2 2a of 1.1 x 10 -3 is 
obtained at Am 2 ~ 300 eV 2 . For sin 2 2a = 1, Am 2 > 1.6 eV 2 is excluded, and for Am 2 > 1000 eV 2 , 
sin 2 2a > 1.8 x 10 -3 is excluded. This result is the most stringent limit to date for Am 2 > 25 eV 2 
and it excludes the high Am 2 oscillation region favoured by the LSND experiment. The u^-to-Ve 
cross-section ratio was measured as a test of <-» universality to be 1.026 ± 0.055. 

PACS numbers: 14.60.Pq, 13.15.+g 
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The existence of neutrino mass and mixing would have 
important implications for fundamental problems in both 
particle physics and cosmology. These include violation 
of lepton family number conservation, the mass of the 
universe, and the observed neutrino deficits from the sun 
and from atmospheric sources. Neutrino oscillations are 
a necessary consequence of non-zero neutrino mass and 
mixing since neutrinos are produced and detected in the 
form of weak-interaction eigenstates whereas their mo- 
tion as they propagate from the point of production to 
their detection is dictated by the mass eigenstates |lj . In 
the two-generation formalism, the mixing probability is: 
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where Am 2 is the mass squared difference of the mass 
eigenstates in eV 2 , a is the mixing angle, E„ is the in- 
coming neutrino energy in GeV, and L is the distance 
between the point of creation and detection in km. 

To date the best limits from accelerator experiments 
for — > i/ e oscillations come from fine-grained calori- 
metric (e.g.: BNL-E734 @, BNL-E776 @) or fully ac- 
tive detectors (e.g. KARMEN [|, LSND §) search- 
ing for quasi-elastic charged current production of elec- 
trons. The LSND experiment, using a liquid scintil- 
lator neutrino target, has reported a signal consistent 
with P M — > D e oscillations at a sin 2 2a w 10~ 2 and 
Am 2 £ 1 eV 2 |§. The CCFR collaboration has previ- 
ously reported a limit on —> i/ e oscillations using the 
ratio of neutral to charged current neutrino events com- 
parable in sensitivity to the above mentioned limits pj . 



In this report we present new limits on ja, 



oscilla- 



tions based on the statistical separation of v e N charged 
current interactions. 

The CCFR detector (7|| consists of an 18 m long, 
690 ton total absorption target calorimeter with a mean 
density of 4.2g/cm 3 , followed by a 10 m long iron toroidal 
spectrometer. The target consists of 168 steel plates, 
each 3m x 3m x 5.15cm, instrumented with liquid scin- 
tillation counters placed every two steel plates and drift 
chambers spaced every four plates. The separation be- 
tween scintillation counters corresponds to 6 radiation 
lengths, and the ratio of electromagnetic to hadronic re- 
sponse of the calorimeter is 1.05. The toroid spectrome- 
ter is not directly used in this analysis which is based on 
the shower profiles in the target-calorimeter. 

The Fermilab Tevatron Quadrupole Triplet neutrino 
beam is a high-intensity, non-sign-selected wideband 
beam with a v : v flux ratio of about 2.5 : 1 and usable 
neutrino energies up to 600 GeV. The production tar- 
get is located 1.4 km upstream of the neutrino detector 
and is followed by a 0.5 km decay region. The resulting 
neutrino energy spectra for f M , v^, v e , and v e induced 
events are shown in Figure [i] The beam contains a 2.3% 
fraction of electron neutrinos, 82% of which are produced 
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FIG. 1. Neutrino energy spectra for v^, i/ e , and v £ at 
the CCFR detector for the FNAL wideband neutrino beam 
(Monte Carlo based on relative and fluxes) . 

from — > n°e ± v e . 

The neutrino interactions observed in the detector can 
be divided into three classes depending on the type of 
the incoming neutrino and on the interaction type: 

1. i/^N — > n~X {Ufj, charged current (CC) events). 

2. v^N — > v^.eX (iSp, e neutral current (NC) events). 

3. v e N -> eX {y e CC events). 

All three types of neutrino interactions initiate a cas- 
cade of hadrons that is registered by the drift chambers 
and scintillation counters. The CC events are char- 
acterized by the presence of a muon produced in the fi- 
nal state which penetrates beyond the end of the hadron 
shower, depositing energy characteristic of a minimum 
ionizing particle |7]] in a large number of consecutive scin- 
tillation counters. Conversely, the electron produced in 
a v e CC event deposits energy in a few counters immedi- 
ately downstream of the interaction vertex which changes 
the energy deposition profile of the shower. The elec- 
tromagnetic shower is typically much shorter than the 
hadron shower and the two cannot be separated for a 
u e CC event. 

In this analysis four experimental quantities are cal- 
culated for each event: the length, the transverse vertex 
position, the visible energy and the shower energy depo- 
sition profile. The event length is determined to be the 
number of scintillation counters spanned from the event 
vertex to to the last counter with a minimum-ionizing 
pulse height. The mean position of the hits in the drift 
chamber immediately downstream of the interaction ver- 
tex determines the transverse vertex position. The visible 
energy in the calorimeter, E V i S is obtained by summing 
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the energy deposited in the scintillation counters from 
the interaction vertex to five counters beyond the end of 
the shower. The shower energy deposition profile is char- 
acterized by the ratio of the sum of the energy deposited 
in the first three scintillation counters to the total visible 
energy. Accordingly, we define 



where Ei is the energy deposited in the i th scintillation 
counter downstream of the interaction place. 

The most downstream counter with energy deposited 
from the products of the neutrino interaction (CEXIT) 
occurs at the end of the hadron shower for NC and 
v e CC events but is determined by the muon track for 
most CC events. We isolate the events without a 
muon track by requiring CEXIT to be no more than 10 
counters downstream from the end of the hadron shower. 
We parametrize the event length which contains 99% of 
such events as: 

L NC = 4. + 3.81 x \og{E ms ) (3) 

In order to measure the number of v e CC events we 
divide the neutrino events into two classes: "short" if 
they deposit energy over an interval shorter than Ljvc, 
and "long" otherwise. The long events consist almost 
exclusively of class 1 events, while the short ones are a 
mixture of class 2, class 3 and class 1 events with a low 
energy muon which cannot be separated on an event-by- 
event basis. 

Based on Lund studies, we take the hadron showers 
produced in NC and CC interactions to be the same. 
Any difference in the shower energy deposition profile 
of long and short events is attributed to the presence 
of v e CC interactions in the short sample. To compare 
directly the long and short events a muon track from the 
data was added to the short events to compensate for the 
absence of a muon in NC events. The fraction, /, of CC 
events with a low energy muon contained in the short 
sample which now have two muon tracks was estimated 
from a detailed Monte Carlo of the experiment in the 
range of 20%. A simulated sample of such events was 
obtained by choosing long events with the appropriate 
energy distribution from the data to which a second short 
muon track was added in software. The length of the 
short track and the angular distribution were obtained 
from a Monte Carlo of CC events. 

To simulate v e interactions in our detector we assume 
Vy, — v e universality. The electron neutrino showers were 
generated by adding a GEANT @ generated electromag- 
netic shower of the appropriate energy to events in the 
long data sample. The energy distribution of the elec- 
tron neutrinos and the fractional energy transfer y were 
generated using a detailed Monte Carlo simulation of the 
experiment. Since the hadron showers in the long sample 




FIG. 2. Eta distributions for short (solid line) , long (dashed 
line) and v e (dotted line) events in four of the energy bins 
studied. The v e and long distributions are normalized to the 
respective number of events predicted by the fit. 

already have a muon track, the v e sample can be com- 
pared directly with the short and long events. 

The long and short % distributions were further cor- 
rected by subtracting the contamination due to cosmic 
ray events. The cosmic ray background was estimated 
from the event sample collected during a beam off gate 
using an identical analysis procedure as for the data 
gates. Additionally, the % distribution of short CC 
events, normalized to the predicted fraction /, was sub- 
tracted from the short event sample. The rjz distributions 
for short, long, and v e events for various energy bins are 
shown in Figure |^. 

For this oscillation search we measure the absolute flux 
of v^s at the detector and compare it to the flux pre- 
dicted by a detailed beamline simulation ||l0| . Any ex- 
cess could be interpreted as a signal of — » v e oscilla- 
tions. The flux was determined directly from the low 
hadron energy CC event sample, normalized to the total 
neutrino cross-section jyl. The same beamline simula- 
tion is used to tag the creation point of each simulated 
Vfi along the decay pipe, and give the number of pre- 
dicted v^s at the detector normalized to the number 
observed at the detector divided by 1 — P(v^ — > v e ). 
P{vn — > v e ) = P(vn —* v e ) is the oscillation probability 
determined from eq. (|l|), assuming CP invariance. The 
predicted electron neutrino flux is normalized to the pro- 
duced number of v^s. The v e flux from neutrino oscil- 
lations is calculated by multiplying the produced number 
Of v^s by P(Uf, -> v e ). 

The events selected are required to deposit a minimum 
of 30 GeV in the target calorimeter to ensure complete 
efficiency of the energy deposition trigger. Additionally, 
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FIG. 3. Number of electron neutrinos as a function of vis- 
ible energy. For electron neutrinos the visible energy is equal 
to the total neutrino energy. The filled band shows Monte 
Carlo prediction assuming no oscillations. The dotted curve 
corresponds to — > v e oscillations with Am 2 = 2000 eV 2 
and sin 2 2a = 0.01 and the dashed curve to Am 2 = 100 eV 2 
and sin 2 2a = 0.01 

we require the event vertex to be more than 5 counters 
from the upstream end of the target and five counters 
plus the separation length from the downstream end and 
less than 50" from the detector centre-line. The resulting 
data sample consists of 632338 long events and 291354 
short ones. 

To extract the number of v e CC events in each of 15 
E V i s bins, we fit the corrected shape of the observed 773 
distribution for the short sample to a combination of 
CC and v e CC distributions with appropriate muon 
additions: 

z^NC(+/z) = a f M CC + /3 z/ c CC(+/i) (4) 

The x 2 °f the fit in each of the 15 E V i S bins ranges from 
33.2 to 77.7 for 41 degrees of freedom (DoF) with a mean 
value of 48.4. Figure || shows that the measured number 
of v e CC's agrees with the Monte Carlo prediction in each 
energy bin. The x 2 value with a no-oscillations assump- 
tion is 9.97/15 DoF. 

The major sources of uncertainties in the comparison 
of the electron flux extracted from the data to that pre- 
dicted by the Monte Carlo are: (i) The statistical error 
from the fit in the extraction of the v e flux, (ii) The error 
in the shower shape modeling, estimated by extracting 
the v e flux using two definitions of ?/. Analogous to the 
definition of 773 given in eq. (|J), we define 774 to be the 
ratio of the sum of the energy deposited outside the first 
four scintillation counters to the total visible energy. If 



TABLE I. The result for sin 2 2a from the fit at each Am 2 
for — > v E oscillations. The 90% C.L. upper limit is equal 
to the best fit sin 2 2a + 1.28a. 



Am (eV ) 


Best fit 


a 


Am (eV ) 


Best fit 


a 


1.0 


-0.1741 


1.6501 


175.0 


0.0000 


0.0016 


2.0 


-0.0501 


0.4107 


200.0 


-0.0002 


0.0014 


3.0 


-0.0153 


0.1852 


225.0 


-0.0003 


0.0013 


4.0 


-0.0112 


0.1041 


250.0 


-0.0004 


0.0012 


5.0 


-0.0051 


0.0671 


275.0 


-0.0004 


0.0012 


7.0 


-0.0036 


0.0345 


300.0 


-0.0004 


0.0012 


9.0 


-0.0021 


0.0213 


350.0 


-0.0004 


0.0012 


10.0 


-0.0023 


0.0173 


400.0 


-0.0003 


0.0013 


20.0 


-0.0004 


0.0048 


450.0 


-0.0003 


0.0015 


30.0 


r\ r\r\r\r\ 

-0.0003 


0.0026 


500.0 


-0.0004 


0.0016 


40.0 


-0.0002 


0.0018 


600.0 


-0.0005 


0.0019 


50.0 


-0.0002 


0.0015 


700.0 


-0.0003 


0.0018 


60.0 


-0.0002 


0.0014 


800.0 


-0.0002 


0.0018 


70.0 


-0.0002 


0.0014 


1000.0 


-0.0004 


0.0017 


80.0 


-0.0003 


0.0014 


1500.0 


-0.0003 


0.0017 


90.0 


-0.0003 


0.0015 


2000.0 


-0.0004 


0.0017 


100.0 


-0.0002 


0.0015 


5000.0 


-0.0003 


0.0018 


125.0 


0.0004 


0.0018 


10000.0 


-0.0004 


0.0017 


150.0 


0.0005 


0.0019 


20000.0 


-0.0004 


0.0017 



the modeling of the showers were correct, the difference 
in the number of electron neutrinos measured by the two 
methods should be small, any difference is used to esti- 
mate the systematic error. Since this error was shown not 
to be correlated among energy bins, we add it in quadra- 
ture to the statistical error from the fit and take this to be 
the combined basic error, (iii) The 1% uncertainty in the 
absolute energy calibration of the detector changes the 
relative neutrino flux which is extracted using the subset 
of the data sample with low hadron energy |[l| . (iv) The 
uncertainty in the incident flux of v^s on the detector is 
estimated to be 4.1% jicj. This error is dominated by 
a 20% production uncertainty in the Kl content of the 
secondary beam which produces 16% of the v e flux. The 
majority of the v e flux comes from i"Trj decays, which are 
well-constrained by the observed spectrum from K^ 2 
decays Jig ]. Other sources of systematic errors were also 
investigated and found to be small. 

The data are fit by forming a x 2 which incorporates the 
Monte Carlo generated effect of oscillations, the basic er- 
ror, and terms with coefficients accounting for systematic 
uncertainties. A best fit sin 2 2a is determined for each 
Am 2 by minimizing the x 2 as a function of sin 2 2a and 
these systematic coefficients. At all Am 2 , the data are 
consistent with no observed — ► v e oscillations. The 
statistical significance of the best-fit oscillation at any 
Am 2 is at most 0.3er. 

The frequentist approach is used to set a 90% con- 
fidence upper limit for each Am 2 . The limit in sin 2 2a 
corresponds to a shift of 1.64 units in x 2 from the mini- 
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FIG. 4. Excluded region of sin 2 2a and Am 2 for — » v e 
oscillations from this analysis at 90% confidence is the area 
to the right of the dark, solid curve. 

mum x 2 (at the best fit value in Table J]) . The 90% confi- 
dence upper limit is plotted in Figure Effor — > v e . The 
best limit of sin 2 2a < 1.1 x 1CT 3 is at Am 2 = 300 eV 2 . 
For sin 2 2a — 1, Am 2 > 1.6 eV 2 is excluded, and for 
Am 2 > 1000 eV 2 , sin 2 2a > 1.8 x 10~ 3 . 

Under the assumption that there are no oscillations, 
this data can also be used to test v^{i>^) <-> v e (p e ) uni- 
versality by comparing the observed v e flux to that pre- 
dicted by the Monte Carlo. From this comparison we 
determine the ratio of the cross sections averaged over 
our flux to be <J C c{ v ii)/ !J cc( l 'e) = 1.026 ±0.055. This is 
currently the most stringent test of universality at high 
space-like momentum transfer. 

In conclusion, we have used the difference in the longi- 
tudinal shower energy deposition pattern of v e N versus 
VpN interactions to search for — > v e oscillations with 
a coarse-grained calorimetric detector. We see a result 
consistent with no neutrino oscillations and find 90% con- 
fidence level excluded regions in the sin 2 2a — Am 2 phase 
space. This result is the most stringent limit to date for 
Vu. ~~ > v e oscillation for Am 2 > 25 eV 2 . We also tested 
^u(^u) v e(pe) universality and found the ratio of the 
v^-io-Vf, cross-section to be 1.026 ± 0.055. 
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